Stability of superfluid and supersolid phases of dipolar bosons in optical lattices 
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We perform a stability analysis of superfluid (SF) and supersolid (SS) phases of polarized dipolar 
bosons in two-dimensional optical lattices at high filling factors and zero temperature, and obtain 
the phase boundaries between SF, checkerboard SS (CSS), striped SS (SSS), and collapse. We show 
that the phase diagram can be explored through the application of an external field and the tuning 
of its direction with respect to the optical lattice plane. In particular, we find a transition between 
the CSS and SSS phases. 
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Since the realization of the superfluid-Mott insulator 
transition of ultracold Bose gases confined to optical lat- 
tices, ultracold atoms have become the playground for 
the realization of various quantum phases studied in con- 
densed matter physics The unprecedented control of 
the lattice depth, dimensionality, geometry, and filling 
factor has allowed for the exploration of a variety of ef- 
fects, including the experimental observation of a dipolar 
condensate of ^^Cr atoms 0,0,01 and the production of 
ultracold heteronuclear molecules 

The problem of interacting dipolar bosons is impor- 
tant not only for magnetic dipolar atoms, but also for 
heteronuclear molecules and Rydberg atoms, which can 
have potentially large electric dipole moments. Thus far, 
dipolar superfluids (SF) have been found only for mag- 
netic dipolar atoms in optical traps 0, 0, Q , but there 
are still several phases that can be pursued experimen- 
tally, including dipolar supersolids (SS), which are char- 
acterized by the simultaneous existence of crystalline and 
SF orders. The possibility of SS phases first emerged 
in the context of solid ^He [^, 0, 0]. More recently, 
there have been experimental reports that the theoret- 
ically predicted 0| non-classical rotational inertia was 
found [9| in solid *He. Although the existence of a SS 
still remains a controversial issue in the condensed mat- 
ter literature [lo| , its existence may be easier to verify in 
the context of Bose gases in optical lattices. 

In this manuscript, we analyze SF and SS phases of 
ultracold dipolar Bose gases loaded into two-dimensional 
(2D) optical lattices, and focus on the region of high fill- 
ing factors. We show that using an external field the sign 
and magnitude of dipole interactions can be controlled 
leading to a variety of diflFerent phases. The phases de- 
scribed include SF, striped SS (SSS), checkerboard SS 
(CSS), and collapse, for which we analyze experimentally 
relevant quantities such as the excitation spectra. 

Our work on ultracold dipolar bosons distinguishes it- 
self from recent work on this topic in several ways. First, 
our discussion is analytical in contrast to numerical work 
based on the Gutzwiller projection techniques_ 11 ^ 12, 13 1 



Second, while some work [H, [13 was confined to purely 
repulsive dipolar interactions, we allow for the compe- 
tition between attractive and repulsive dipolar interac- 
tions by changing the direction of an external field with 
respect to the 2D lattice plane. Third, we explore the 
phase diagram for a wider range of interactions than pre- 
viously investigated 



111, liJ, [13|, [17| and find that the 



quantum phase of the system can be switched between 
CSS and SSS as the direction of the dipole is changed. 
The transition between these SS phases is remarkable in 
the sense that it is a structural phase transition between 
two types of SS, where the symmetry of the crystalline 
order changes. Lastly, we go beyond mean-field descrip- 
tions by including fluctuation effects and performing a 
stability analysis. 

To study bosons with dipole-dipole interactions in 2D 
optical lattices, we use the dipolar-Bose-Hubbard model 



H^-jJ2{b]bi+h.c.) + -J2n,{n, 



-1) 



^Y^Vjinjni, (1) 

3<l 



where b 



is the boson creation operator at site j, 
and J is the hopping. (j7) represents 



nearest-neighbor (NN) pairs of lattice sites. The on- 
site interaction U = Us + Ud consists of the con- 
tributions from both the s-wave scattering Us = 
47rS^as/r7i J dr|?i'(r)|^ and the onsite dipole-dipole inter- 



or quantum Monte Carlo (QMC) methods [llTlSl. 11^. 



action Ud = (27r)3/2/ dkn^{k)Vd{k) where is 

the s-wave scattering length, w{r) is the Wannier func- 
tion with respect to the underlying lattice potential, and 
n(k) and Vd(k) are the Fourier transform of the density 
|w(r)p and the dipole-dipole potential. The long-range 
part of the dipole-dipole interaction is well-approximated 
as Vji = D^{l~3cos^9ji)\rj — ri\~^ where D is the dipole 
moment and Oji is the angle between the dipole direc- 
tion and Tj — ri. Here, = {j^a^jyO) is a lattice vector, 
where jx and jy are integers and a is the lattice spacing. 
A schematic picture of the system is shown in Fig. [TJ 

Since we are interested only in SF and SS phases, we 
restrict ourselves to Bose-condensed solutions of Eq. ([IJ . 
Minimizing the quantum action S = ih^^ b^jdtbj — H 



2 



with respect to the saddle-point field j (t) obtained from 
the transformation bj{t) ^jit) = ^j{t) + 
and neglecting the fluctuations <5^j(t) lead to the time- 
dependent Gross-Pitaevskii (GP) equation ihdt'^fj = 

-JEn) + {U\^A^ + Ei^j V,i\^i\')^,. The dcscrip- 
tion above is justified when the healing length ^ = 
h/(m*c) is larger than the interatomic distance / = 
aj^Jv [l9|, where c is the sound velocity, m* is the ef- 
fective mass, and v is the filling factor. This condition 
is satisfied when J 3> max([/, In addition, as dis- 

cussed later, the emergence of SS phases requires J < v\J ^ 
which can be satisfied when i/ ^ 1. Thus, our approach 
is valid for large hopping and large filling factor. Notice 
that we are not investigating incompressible solid phases, 
such as Mott-insulators and density wave solids, because 
they can not be described within this GP approximation. 

Writing ^'j(t) = $^e^'^*/'' requires the static part of 
the condensate to satisfy 

- J^ci>, + (C/|<f>,f + = M*,, (2) 

where /i is the chemical potential. The normalization 
condition is M~^^^ I'J'jP = where M is the total 
number of sites; and the energy of the condensate is 

m 3 3<i 



A standard procedure to obtain phase diagrams is 
the comparison of energies for various phases as done 
in mean-field studies based on the Gutzwillcr projec- 
tion technique [ll|, [3]. It is also important to 
study the stability of each phase in order to construct 
the correct phase diagram. For this purpose, we perform 
a stability analysis by taking into account fluctuations 
<5*j- = e-^'^*/'*(uje-*^*/'* - v*e''^'^/^) beyond the saddle- 




FIG. 1: (color online) (a) Schematic picture of a 2D lattice 
system, where a square represents a site, (b) The dipole vector 
is shown in the spherical polar coordinates. 



point ^1 j, which lead to the eigenvalue equations 

(/> ^ 1^3 ^ 

E '^3,i^*i^i-U^']v, - $j E V,- =eu,, (4) 
1^3 17^3 

-jJ2^l+(2U\<i>3\' + T.^3A^'\^-Av, 

(0 ^ 1^3 ^ 

1=^3 1^3 

Here, e is the energy and {uj,Vj) is the amplitude of 
the collective mode of the condensate. The appearance 
of collective modes with complex frequencies signals ex- 
ponential growth of fluctuations in time, and thus the 
existence of dynamical instabilities (DI). 

Since the dipolar interaction decays over distance, 
the ratios between the NN interactions I4 = (1 — 
3sin^acos2 0)D2/a3, Vy = (1 - 3 sin^ a sin^ 0)1)2/0^ 
(along the x and y directions) and the onsite interac- 
tion U dictates the basic physics. Here, a and are the 
elevation and azimuthal angles of the direction of the po- 
larization with respect to the 2D lattice plane as shown in 
Fig. Hl^b). These ratios are experimentally controllable, 
e.g., by changing the s-wave scattering length via a Fes- 
hbach resonance, as demonstrated with ^^Cr atoms 
Moreover, the ratio 7 = Vy/Vx is also controllable by 
changing the polarization direction through the applica- 
tion of an external fleld. For instance, when (j) = tt/2 and 
a varies, the relevant ratio becomes 7 = 1 — Bsin^a. 

To gain analytical insight, we consider first onsite and 
NN interactions. In this case, there are four possible 
pure condensate phases. When the absolute values of Vx 
and Vy are much smaller than U, the system is in a SF 
phase. When the NN interactions are almost isotropic 
and strongly repulsive, diagonal crystalline order devel- 
ops and a CSS phase emerges. When the NN interaction 
is strongly anisotropic, e.g. Vx = —Vy, a SSS phase is 
favored. Schematic pictures of these phases are shown in 
Figs. [2i;i)-(IV). The phase boundaries of SF, CSS, and 
SSS can be determined by solving Eq. ([2]) and obtaining 
the ground state energy given in Eq. ([3|). In addition, 
strongly attractive NN interactions result in the collapse 
of the condensate, which is characterized by a DI occur- 
ring in the long-wavelength phonon mode. In Fig. [21 the 
solid and dotted lines indicate the phase boundaries in 
the {Vx/U, Vy/U) plane for J/{iyU) = 0.1. 

Let us first identify the SF region. In the SF phase, the 
condensate wave function, $j = is uniform and the 
chemical potential is /isF = —4:J+i^{U +2Vx+2Vy) . From 
Eq. Q, we obtain the energy Esf/M = -AJv+v'^{U /2+ 
Vx + Vy), and from Eqs. (j4|) and ([5]) we obtain the excita- 
tion spectrum esF(q) = \/e(q) [e(q) + 2i/(t/ + 2F(q))], 
where q is the quasimomentum of the collective mode. 
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FIG. 2: (color online) Left: Phase diagrams in the 
(y:^/U,Vy/U) plane for J/{vU) = 0.1. While only the NN 
interactions are considered in (a) , the full dipolar interactions 
are included in (b). The solid (dotted) lines represent phase 
boundaries for continuous (discontinuous) transitions. The 
dashed-dotted green lines and the dashed red lines separate 
the stable and unstable regions of the CSS state and those of 
the SSS. Right: Schematic pictures of SF (I), CSS (II), SSSl 
(III), SSS2 (IV), SSal (V), SSbl (VI), and SScl (VII). 



g(q) = 4J (sin^(g:ra/2) + sin^((7j^a/2)), and F(q) = 
VxCOs{qxa) + VyCos{qya). 

Since esF(q) — csFfiq for qa <^ 1, where csf = 
\j2viJJ + 2Vx + 2Vy)Ja/h is the sound velocity for the 
SF phase, long- wavelength phonons cause a DI when 
U + 2Vx + 2Vy < for any non-zero value of J. The sound 
velocity c is directly related to the compressibility k via 
= m*c^, which clearly shows that the compressibil- 
ity becomes negative at this DI leading to the collapse of 
the condensate [2^. On the other hand, since EsfCQo) = 
^/8IT2^iU^^Wx^^2yy), where Qo = (7r/a,7r/a), the 
collective modes in the vicinity of Qo cause a DI when 
214 -I- 2Vy — U — AJ/v > 0, signaling a transition to the 
CSS phase [l3|, where the mode with q = Qo is am- 
plified and its interference with the condensate in the 
zero momentum state creates the checkerboard density 
wave order. Similarly, we find that the collective modes 
in the vicinity of q = Qi = (tt/o, 0) or Q2 = (0,7r/a) 
cause a DI and signal transitions to the SSS phases, when 
2Vx -2Vy-U- 2J/i> > or 2Vy - 2Vx ~U - 2J/v > 0, 
respectively. Thus, the SF phase corresponds to the 
square region surrounded by the lines 214 + 2Vy + U ^ 0, 
2T4 + 2Vy-U - AJ/v = 0, 214 - 2Vy-U - 2J/v = 




FIG. 3: (color online) Excitation spectra ecgg(ga:, % — 0) and 
^sssi(9^ = 0,qy) for J/{uU) — 0.1, and 14 — U. (a) The real 
part and (b) the imaginary parts of e^gg {qx ,qy = 0) . (c) The 
real part and (d) the imaginary parts of eggsi(?a; = 0, Qy). 

and 2Vy - 2^ - [/ - 2 J/j/ = as shown in Fig. E^a). 

Since the phase diagram is symmetric with respect 
to the line 14 = Vy, we focus on the region Vx > Vy 
henceforth. The condensate wave functions for the CSS 
and SSSl states are <^f^^ = ^/i^ + yr^e**^" '"^ and 
•fj^^^ = y/t^ + y/i^e''^^'''^ . Here, z^c {i^s) is the den- 
sity of atoms condensed in the state with quasimomen- 
tum Qo (Qi), corresponding to the order parameter of 
CSS (SSSl). These states are regarded as SS phases in 
the sense that each of them possesses off-diagonal long- 
range order through the condensate, and diagonal long- 
range order via a density wave. From Eq. for the 
CSS phase we obtain ^css = 2i/C/, i^co = v ~ v^, and 
i/c = yl2 - 2J/[2Vx + 2Vy - U). Similarly for SSSl, we 
obtain ^sssi = —2J + 2v{U + 2Vy), Vso = v — Vg, and 
i^s ~ iy/2 — J/{2Vx — 2Vy — U). Since Vc {vs) vanishes 
at the phase boundary between SF and CSS (SSSl), the 
phase transitions are continuous according to Landau's 
classification (or second-order according to Ehrenfcst's). 

From Eq. ([3]), the energies of the condensate for 
the CSS and the SSSl are expressed as i?cssA^ = 
-8JV(214 + 2Vy -U) + ly^U, and £;sssi/^/ = -2Jiy - 
2JV(2K - 2Vy - U) + v'^{2Vy + U). The condition, 
-E-css = -E-sssij determines the phase boundary between 
CSS and SSSl as shown in Fig. (S^a). Both in CSS 
and Us in SSSl are finite at the phase boundary. This 
fact indicates that the phase transition between CSS 
and SSSl is discontinuous (first-order). A discontinu- 
ous phase transition is characterized by hysteresis, such 
that the critical point Vy^ of the transition from CSS to 
SSSl disagrees with the critical point Vy'^ of transition 
from SSSl to CSS due to the presence of more than one 
meta-stable state. To determine Vy^ and Vy'^ , we calcu- 
late the excitation spectra ecss(^) ^'^'^ ^sssi 

(q) for CSS 

and SSSl, where the plus (minus) sign corresponds to a 
gapful (gaplcss) mode at low momenta. 
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In Figs. [Slja) and (b), we show £cssilx^1v ~ ^) ^'^^ 
varying Vy, but fixed = U. As decreases, a roton- 
like minimum is formed at q = Qi and reaches zero at 
Vy = Vy^ . As is decreased further, the imaginary part 
of ecss(Qi) grows, thus revealing the transition to the 
SSSl phase. The condition ecss(Qi) = ^ ^^"^'^ gives Vy^ 
as shown by the dashed-dotted green hue in Fig.[2Ua). In 
Figs. [3^c) and (d), we show £ggsi('?2; = 0,%) for varying 
Vy, but fixed Vx = U. The condition £gggi(Q2) = gives 
Vy^ as shown by the dashed red Hne in Fig. [Slja). 

To complete our analysis of the phase diagram, we lo- 
cate now the boundary between phases SSSl and col- 
lapse. When qa 1, 



esssM -V^Vy + U 



2K 



-2Vy-U 



-4:iyJ{qy 



1/2 



,(6) 



and the DI leads to the collapse when 2Vy + U < (sec 
Fig. Ha)). 

Next we numerically calculate the energy and per- 
form linear stability analyses for the full dipolar in- 
teractions. We choose a specific case of ~ or 
7r/2, where the dipolar interactions can be written as 
Vji = {ilVx + ilVy)/{il + ix and iy are inte- 

gers that satisfy Vj — r; = {ixa,iya). In this case, we 
obtain the phase diagram in the (T4/C/, Vy/i7)-plane for 
J/^uU) = 0.1 as shown in Fig.[2lb) [2]J. As well as in the 
phase diagram with only the NN interactions, there re- 
main the CSS and SSS regions, which we have confirmed 
to be stable. CSS shares phase boundaries with SSS and 
the transition between these two phases is discontinuous 
(first-order). However, unlike the case with only the NN 
interactions, there are additional SS phases with differ- 
ent density wave orders, such as SSa, SSb, and SSc (see 
(V)-(VII) of Fig. [2). While SSb has two-sublattice den- 
sity modulation as well as CSS and SSS do, SSa and SSc 
have three-sublattice modulation. There can be other 
SS phases with more than three-sublattice modulation, 
which shares boundary with SSa and SSc, but here we 
do not determine the location of these phases because we 
are mainly interested in the SF. CSS, and SSS phases. 

In phenomenological classical models involving dipo- 



24| 



lar interactions for 2D continuum systems |22l . [2 
it has been conjectured that superstructures or micro- 
emulsion phases exist between two stable phases in a re- 
gion of the temperature versus density phase diagram 
where phase separation overrides an otherwise expected 
first order phase transition. A necessary condition for the 
emergence of such micro-emulsions at zero temperature 
is the negative compressibility leading to phase separa- 
tion. For a first order transition between CSS and SSS 
described within our model and approach, we do not find 
a negative compressibility, namely the pre-condition for 
phase separation, in the transition region. In contrast, 
recent QMC simulations by PoUet et al. have shown pre- 



liminary evidence of micro-emulsion phases in systems of 
hardcore dipolar bosons in a 2D triangular lattice [1^. 
However, QMC simulations of larger system sizes are nec- 
essary to make a conclusive statement about the presence 
of the micro-emulsion phases. 

In conclusion, we have studied stability of quantum 
phases of dipolar bosons in two-dimensional optical lat- 
tices. We have shown that mean-field theories alone fail 
to describe correctly the phase boundaries, and that a 
stability analysis is absolutely necessary. In addition 
to superfiuid and collapsed phases, we have shown that 
striped and checkerboard supersolids can exist, and com- 
pete with each other due to the anisotropy of dipolar 
interactions as controlled by an external field. 
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thanks NSF (DMR-0709584) and ARO (W911NF-09-1- 
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